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ABSTRACT 

Photometric and spectroscopic data of the energetic Type Ic supernova (SN) 2002ap are 
presented, and the properties of the SN are investigated through models of its spectral 
evolution and its light curve. The SN is spectroscopically similar to the "hypernova" 
SN 1997ef. However, its kinetic energy [~ (4 — 10) x 10^^ erg] and the mass ejected (2.5-5 
Mq) are smaller, resulting in a faster-evolving light curve. The SN synthesized ~ O.OTMq 
of ^^Ni, and its peak luminosity was similar to that of normal SNe. Brightness alone 
should not be used to define a hypernova, whose defining character, namely very broad 
spectral features, is the result of a high kinetic energy. The likely main-sequence mass of 
the progenitor star was 20-25 Mq, which is also lower than that of both hypernovae SNe 
1997ef and 1998bw. SN 2002ap appears to lie at the low-energy and low-mass end of the 
hypernova sequence as it is known so far. Observations of the nebular spectrum, which 
is expected to dominate by summer 2002, are necessary to confirm these values. 

Subject headings: supernovae: general - supernovae: individual (SN 2002ap) - line: 

identification - line: formation - nucleosynthesis - gamma rays: bursts 
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1. INTRODUCTION 

One of the most interesting recent developments 
in the study of supernovae (SNe) is the discov- 
ery of some very energetic Type Ic SNe (SNe Ic; 
see Fihppenko 1997 for a general review), whose 
kinetic energy (KE) exceeds 10^^ erg, about 10 
times the KE of normal core-collapse SNe (here- 
after E^i = 10^^ erg). The most luminous and 
powerful of these objects, SN 1998bw, was prob- 
ably linked to GRB 980425 (Galama et al. 1998), 
thus establishing for the first time a connection 
between the enigmatic gamma-ray bursts (GRBs) 
and the well-studied phenomenon of core-collapse 
SNe. However, SN 1998bw was exceptional for a 
SN Ic: it was as luminous at peak as a SN la, in- 
dicating that it synthesized ~ 0.5 Mq of ^^Ni, and 
its KE was estimated at ~ 3 x 10^^ erg (Iwamoto 
et al. 1998; Woosley et al. 1999). Because of its 
large KE, SN 1998bw was classified as a "hyper- 
nova". Signatures of asymmetry were detected 
in SN 1998bw, such as polarization (Patat et al. 

2001) and pecuhar nebular-line profiles (Mazzali 
et al. 2001; Macda et al. 2002). 

Subsequently, other "hypernovae" of Type Ic 
have been discovered or recognised, such as 
SN 1997ef (Iwamoto et al. 2000; Mazzah et al. 
2000) and SN 1997dq (Matheson et al. 2001; Maz- 
zali et al. 2002, in preparation), although their 
KE was not as large as that of SN 1998bw, 
and they did not appear to be associated with 
GRBs. The analysis of these various objects sug- 
gests that the KE may be related to the progeni- 
tor's main-sequence mass, which was ~ 40 Mq for 
SN 1998bw and ~ 30 Mq for SN 1997ef. These 
values place hypernovae at the high-mass end of 
SN progenitors. Another possible hypernova, al- 
though of Type Iln, was SN 1997cy, which was 
also estimated to have a large mass (~ 25 Mq; 
Germany et al. 2000; Turatto et al. 2000). It is 
not yet clear whether the large mass is the dis- 
criminating factor for the birth of a hypernova or 
the connection with a GRB. 

In this paper we present a first analysis of the 
properties of the recently discovered SN Ic 2002ap 
in M74 (Hirose 2002). The SN was immediately 
recognised hypernova from its broad spectral 
features (Kinugasa et al. 2002; Meikle et al. 2002; 
Gal- Yam et al. 2002a; Filippenko & Chornock 

2002) . This indicates high velocities in the ejected 
material, which is the typical signature of hy- 
pernovae. It was therefore followed from several 



observatories, and the relative proximity also fa- 
vored observations with small telescopes. Luck- 
ily, the SN was discovered very soon after it ex- 
ploded: the discovery date was January 29, while 
the SN was not detected on January 25 (Nakano 
et al. 2002). This is among the earliest any SN 
has been observed, with the obvious exceptions of 
SN 1987A and SN 1993J. 

Figure 1 shows the near-maximum spectra of 
SN 2002ap, of the hypernovae SNe 1998bw and 
1997ef, and of the normal SN Ic 19941. If hue 
width is the distinguishing feature of a hypernova, 
then clearly SN 2002ap is a hypernova, as its spec- 
trum resembles that of SN 1997cf much more than 
that of SN 19941. Line blending in SN 2002ap and 
SN 1997ef is comparable. However, some individ- 
ual features that are clearly visible in SN 19941 
but completely blended in SN 1997ef can at least 
be discerned in SN 2002ap (e.g. the Na I-Si II 
blend near 6000 A and the Fe ll lines near 5000 A). 
Therefore, spcctroscopically SN 2002ap appears 
to be located just below SN 1997ef in a "velocity 
scale," but far above SN 19941. 

This appears to be confirmed by the light curve. 
Figure 2 shows the F-band light curves of the 
same four SNe as in Figure 1. SN 2002ap reached 
V maximum on about February 8 at V = 12.3 
mag. SN 2002ap peaks earlier than both hyper- 
novae 1998bw and 1997ef, but later than the nor- 
mal SN 19941, suggesting an intermediate value 
of the ejecta mass Mcj. 

Using a distance to M74 of 8 Mpc (// = 29.5 
mag; Sharina et al. 1996), and a combined Galaxy 
and M74 reddening of E{B — V) — 0.09 mag (es- 
timated from a Subaru HDS spectrum; Takada- 
Hidai et al. 2002, in preparation), the absolute 
magnitude is My = —17.4. This is comparable 
to SN 1997ef and fainter than SN 1998bw by al- 
most 2 mag. Since peak brightness depends on the 
ejected ^^Ni mass, SNe 2002ap, 1997ef, and 19941 
appear to have synthesized similar amounts of it. 
Estimates were ~ 0.07 Mq for SN 19941 (Nomoto 
et al. 1994) and 0.13 Mq for SN 1997ef (Mazzah 
et al. 2000). 

Section 2 of this paper describes how an ex- 
plosion model was selected that gives reasonable 
fits to the spectra of SN 2002ap, thus establishing 
the properties of the SN. In § 3 we show synthetic 
light curves obtained with our best models, which 
confirm the results of the spectroscopic analysis. 
In § 4 we discuss implications for the progenitor 
of SN 2002ap. 



3 



2. SPECTROSCOPIC MODELS 

Iwamoto et al. (1998) showed that synthetic 
hght curves of SNe Ic are degenerate with respect 
to particular sets of the parameters Mcj, KE^ and 
the opacity. Spectra must also be used to de- 
termine those parameters, as their dependence is 
different. Therefore we started the analysis from 
the spectra. We used the Monte-Carlo SN spec- 
trum synthesis code described by Mazzali & Lucy 
(1993) and modified as in Lucy (1999) and Maz- 
zah (2000). 

We tested several versions of the hypcrnova 
explosion models CO 138 (used for SN 1998bw; 
Iwamoto et al. 1998) and COlOO (used for 
SN 1997ef; Mazzah et al. 2000). Both of the 
models turned out to be too massive. In fact, 
since SN 2002ap rises very rapidly in brightness, 
the epochs of its spectra are smaller than those 
of similar-looking spectra of both SNe 1997ef and 
1998bw. In particular, the earliest spectrum of 
SN 2002ap, taken on January 30, has an epoch 
of at most 4 days, but it resembles the spectra 
of SNe 1997ef and 1998bw at epochs of about 8 
days, indicating similar temperatures and expan- 
sion velocities. The densities in the massive mod- 
els at such an early epoch are much too high, and 
no good solution can be obtained unless most of 
the mass is composed of He, which contributes 
very little both to the electron density and to line 
opacity. 

As the subsequent evolution of SN 2002ap more 
closely resembles that of SN 1997ef (albeit at a 
faster rate) than that of SN 1998bw, we selected 
model COlOO and rescaled it to lower masses, 
and hence lower KE. We tried a model with 
half the mass (model COlOO/2: Mej= 4.8 Mq; 
i?5i = 8) and one with a quarter of the mass 
(model COlOO/4: Mej= 2.4 M©; E^i = 4). Model 
COlOO/2 is too massive, while model COlOO/4 
appears to be roughly appropriate for SN 2002ap. 
Figure 3 shows a series of fits obtained with that 
model. 

Also, the earliest spectrum (Meikle et al. 2002) 
requires a large degree of line blending, as shown 
by the apparent lack of any separation between 
O I A7774 and the Ca ll near-infrared triplet. 
This requires the presence of sufficient material at 
V > 30, 000 km s~^ In the case of SN 1997ef, this 
blending was obtained by adding an outer "flat" 
part to the density profile (Mazzali et al. 2000). 
Blending in the O l-Ca ll feature in the earliest 



observed spectrum of SN 1997ef line was weaker, 
but the epoch of this spectrum is later than the 
first spectrum of SN 2002ap. We also modified 
model COlOO slightly by limiting the velocities 
to 65,000 km s~^, to match the observations. 

Because the first spectrum was obtained so 
early, it was possible to establish its age fairly 
accurately using synthetic spectra. Age and ve- 
locity must in fact combine to give an acceptable 
temperature, but velocity can be checked against 
the blueshift of the line absorptions. For January 
30 we find t = 2.Q ± 0.5 days. 

The first spectrum has a photospheric velocity 
Vph — 30, 000 km s~^, and thus the photosphere is 
located just below the outer flat extension of the 
density proflle. Hence line blending is very strong, 
and the spectrum resembles that of SN 1998bw. 
But already the next spectrum, with t = 3.5 days, 
has fph = 20,500 km s^^, which is in the steep 
part of the density distribution, so blending be- 
comes less severe. Note that the Ca ll near-IR 
triplet can blend with the O I line even with a very 
small abundance of Ca at high velocity, because it 
is a very easy line to excite. We see that Vph con- 
tinues to decrease with time, but the rate of de- 
crease becomes smaller as the light curve reaches 
peak, which is expected as the photosphere then 
enters the inner flat part of the density distribu- 
tion. 

Oxygen dominates the composition in all the 
spectra we have modeled. However, the earliest 
spectrum contains ~ 30% He by mass, suggesting 
that the SN exploded with about 0.05 M© left 
from the He envelope. This appears to be con- 
sistent with the early detection of an absorption 
at about 1.03 /im, which could be He I A10830, 
although Mg ii and C I could also contribute to it 
(Fig. 3; Motohara et al. 2002). Synthetic spectra 
can only reproduce this feature if signiflcant non- 
thermal ionization of He is active, which is also 
not unexpected. The abundances of Si and other 
intermediate-mass elements, and those of the Fe- 
group elements, increase with depth in the ejecta, 
which is in accord with expectations. 

Several problems remain in the models. The 
synthetic O I and Ca II absorptions are too blue 
at the more advanced epochs, suggesting that the 
density profile we have adopted may overestimate 
the mass at high velocities. Also, the ratio of the 
Fe II lines in the absorption near 4800 A is not 
correct. Further, more refined analysis is there- 
fore needed. However, we do not expect that our 
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general results will change significantly. 

3. LIGHT-CURVE MODELS 

We constructed the UVOIR bolometric 
light curve of SN 2002ap mainly from the 
UBV{RI)cJHK photometry obtained with the 
MAGNUM telescope (Yoshii et al. 2002). We 
also included 5 data points based on observations 
from Wise Observatory (Gal- Yam et al. 2002b), 
the U.K. Infrared Telescope (Mattila et al. 2002), 
MDM (Gerardy et al. 2002, in preparation) and 
the Subaru Telescope (Motohara et al. 2002). We 
estimate the errors at 0.3 mag, mostly owing to 
the uncertainty in the distance modulus (~ 0.2 
mag) and the nonstellar nature of the spectrum. 
The UVOIR magnitudes and the bolometric mag- 
nitudes derived from the synthetic spectra agree 
to within 0.1 mag. The IR flux accounts for a 
significant fraction of the total flux, ranging from 
~ 20% on February 4 to ~ 40% on March 10, 
while the UV flux observed by XMM-Newton on 
February 3 (Pascual et al. 2002) contributes only 
-4%. 

Synthetic light curves were computed with an 
LTE radiation hydrodynamics code and a gray 7- 
ray transfer code (Iwamoto et al. 2000). TOPS 
opacities (Magee et al. 1995) were fitted to the 
ejecta models in order to find an empirical re- 
lationship between the Rosseland mean and the 
electron-scattering opacity. 

Figure 4 shows the bolometric data and two 
synthetic light curves computed with the same 
density structures used for spectrum synthesis. 
The two models produce very similar light curves, 
demonstrating the parameter degeneracy (§2). 
Both model fight curves give a satisfactory fit 
to the observations, confirming the results of the 
spectral analysis. 

To reproduce the observed peak luminosity, a 
total of 0.07 M0of ^^Ni is used. In order to achieve 
a rapid rise of the light curve, and in particu- 
lar to reproduce the earliest bolometric point on 
day 2, we had to mix ^^Ni out to high velocities, 
which introduces a local concentration of ^^Ni 
very close to the surface, as in the best-fit model 
for SN 1998bw (Nakamura et al. 2001). Such sig- 
nificant outward mixing of the newly synthesized 
^^Ni is consistent with our spectral analysis. 

Considering the uncertainties in the UVOIR 
bolometric light curve, we suggest that M (^^Ni) 
= 0.07 ± 0.02 M0. One should not be misled by 



the fact that SN 2002ap is more luminous in the V 
band than SN 19941 to conclude that the former 
has a larger ^^Ni mass. In fact, the bolometric 
corrections (B.C. = Mboi - My) in SN 2002ap 
are quite large, ranging from 0.5 to 0.3, during 
the first 20 days. In contrast, SN 19941 shows 
a small B.C. of 0.05 at the luminosity maximum 
(Iwamoto ct al. 1994). Large values for the B.C. 
are also found in SN 1997ef (Mazzafi et al. 2000). 

4. DISCUSSION AND CONCLUSIONS 

The spectral evolution of SN 2002ap appears to 
follow closely that of SN 1997ef, at about twice 
the rate. The spectra and the light curve of 
SN 2002ap can be well reproduced by a model 
with ejected heavy-element mass Mej= 2.5-5 Mq 
and £^51 = 4-10. Both M^j and KE are much 
smaller than those of SNe 1998bw and 1997ef (but 
they could be larger if a significant amount of He 
is present). The ^^Ni mass is estimated to be 
~ 0.07 M0, which is similar to that of normal 
core-collapse SNe such as SNe 1987A and 19941. 

Although SN 2002ap appears to lie between 
classical core-collapse SNe and hypernovae, it 
should be regarded as a hypernova because its 
kinetic energy is distinctly higher than for clas- 
sical core-collapse SNe. In other words, the 
broad spectral features that characterize hyper- 
novae are the results of a high kinetic energy. 
Also, SN 2002ap was not more luminous than nor- 
mal core-collapse SNe. Therefore brightness alone 
should not be used to discriminate hypernovae 
from normal SNe, while the criterion should be 
a high kinetic energy. Further examples of hyper- 
novae are necessary in order to establish whether 
a firm boundary between the two groups exists. 

For these values of KE, Mq, and M{^^m), 
we can constrain the progenitor's main-sequence 
mass and the remnant mass Mj-em- Modeling 
the explosions of C-l-0 stars with various masses, 
we obtain M(^^Ni) function of the parameter 
set {KE, Mco, M^em =Mej-Mco). The model 
which is most consistent with our estimates of 
(Mcj, KE) is one with Mco ~ 5 My. Mrcm ~ 
2.5 Mq, and ^51 = 4. The 5.0 Mq C+0 core 
forms in a He core of mass M^ = 7.0 Mq, corre- 
sponding to a main-sequence mass M^g ~ 20-25 
M0. The — Ma relation depends on convec- 
tion and metallicity (e.g., Nomoto & Hashimoto 
1988; Umeda & Nomoto 2002). Given the non- 
detection of the progenitor in pre-discovery im- 
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ages of M74 (Smartt et al. 2002), the above mass 
range seems to be most consistent with the pro- 
genitor being a member of an interacting binary 
system. 

In order to determine the mass of O and that 
of ^^Ni, whose decay to ^^Fe via ^^Co will be pow- 
ering the nebular spectrum, it is very important 
that spectra arc taken in the fully nebular epoch. 
There is already some evidence for net emission 
in the Ca ll near-IR triplet in a February 21 spec- 
trum, which has an epoch of 24 days and looks 
similar to SN 1997ef on about day 42. Eventually, 
though, O I A6300 should become the strongest 
line. Because of the fast evolution of SN 2002ap, 
the time when the SN is again visible (late June 
2002), corresponding to an epoch of 5 months, 
should already be suitable. 

The estimated progenitor mass and explosion 
energy are both smaller than those of previous 
"hypernovae" such as SNe 1998bw and 1997ef, 
but larger than those of normal core-collapse SNe 
such as SN 1999em. There appears to be a cor- 
relation between initial stellar mass and explo- 
sion energy, but further observational examples 
are needed to estabhsh the relation. In partic- 
ular, it is unclear what fraction of massive stars 
with M > 20 M0 explode energetically. Massive 
core-collapse SNe with either a normal explosion 
energy (e.g., SN 1984L; Swartz & Wheeler 1991) 
or a very small one (SN 1997D; Turatto et al. 
1999) also appear to exist. 

Given the estimated mass of the progenitor, bi- 
nary interaction including the spiral-in of a com- 
panion star (Nomoto et al. 2001) is probably re- 
quired in order for it to lose its hydrogen and some 
(or most) of its helium envelope. This would sug- 
gest that the progenitor was in a state of high 
rotation. It is possible that a high rotation rate 



and/or envelope ejection are also necessary con- 
ditions for the birth of a hypernova. 

The high explosion energy of SNe 1998bw and 
1997ef may be related to the formation of a black 
hole through the extraction of energy from the 
black hole-accretion disk system. Also in the case 
of SN 2002ap, the expected remnant mass (~ 2.5 
M0) exceeds the maximum neiitron star mass, 
even though the estimated main-sequence mass 
(20-25 Mq) is not extremely large. 

SN 2002ap was not apparently associated with 
a GRB. This may actually be not so surprising, 
since the explosion energy of SN 2002ap is about a 
factor of 5-10 smaller than that of SN 1998bw, as 
also indicated by the weak radio signature (Berger 
et al. 2002). The present data show no clear signa- 
ture of asymmetry, except perhaps for some polar- 
ization (Wang et al. 2002; Kawabata et al. 2002), 
which is smaller than that of SN 1998bw. This 
suggests that the degree of asphericity is smaller 
in SN 2002ap and that the possible "jet" may have 
been weaker, which makes GRB generation more 
difficult. It will therefore be interesting to see 
whether the photosphere persists to very low ve- 
locities, as it did in SN 1997ef, in apparent contra- 
diction with one-dimensional explosion models. 
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Fig. 1. — The near-maximum spectra of Type Ic SNe and hypernovae. 
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Fig. 2. — The observed y-band light curves of SNe 2002ap (stars), 1998bw {open circles), 1997ef {open triangles), and 
19941 {filled circles). Most of the SN 2002ap data are taken from the MAGNUM telescope (Yoshii et al. 2002) and the 
rest from the VSNET ( http : / / www. kusastro . kyoto-u .ac.jp/ vsnct /[ ) . 
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Fig. 3. — A comparison between some observed spectra of SN 2002ap (thick lines: January 30 WHT (Meikle et 
al. 2002); February 1 — Gunma Obs. (Kinugasa et al. 2002); February 6 — Beijing Obs. (Qiu et al., in preparation); 
February 8 — Subaru FOG AS (Kawabata et al. 2002); February 7 IR — Subaru GISGO (Motohara et al. 2002), shown 
twice; February 10 — Lick Obs. (Filippenko et al. 2002); February 21 — Asiago Obs. (Turatto et al., in preparation); ) 
and synthetic spectra computed with model GOlOO/4 {dashed lines). 
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Fig. 4. — The bolometric light curve of SN 2002ap constructed using mainly MAGNUM UBV{RI)cJHK photometry 
{filled circles] five points are based on photometry at Wise Observatory (Gal- Yam et al. 2002b), the U.K. Infrared 
Telescope (Mattila et al. 2002), MDM (Gerardy et al., in preparation) and Subaru Telescope (Motohara et al. 2002)). 
This is compared with two best-fit model light curves (Mej = 2.5 M© and £^51 = 4: red line; Moj = 5 Mq and £^51 = 8: 
green line) . The open squares are luminosities from the synthetic spectra. 



